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Abstract 

We study the consequences of the modification of the phase space structure of General Rela- 
tivity imposed by breaking the full diffeomorphism invariance but retaining the time foliation 
preserving diffeomorphisms. We examine the different sectors in phase space that satisfy the 
new structure of constraints. For some sectors we find an infinite tower of constraints. In 
spite of that, we also show that these sectors allow for solutions, among them some well 
known families of black hole and cosmologies which fulfill all the constraints. We raise some 
physical concerns on the consequences of an absolute Galilean time, on the thermodynamical 
pathologies of such models and on their unusual vacuum structure. 
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1 Motivation, outlook and conclusions 



There has been recently a considerably huge activity regarding different aspects of an even- 
tually renormalizable gravitational theory [1]. Despite the immense amount of work that 
the subject has triggered, only a few authors have dealt in detail with the consistency of 
the initial proposal [21 El IH while the vast majority of them deal with applications in 
cosmology [6l [7] and with the obtention of some solutions, typically applicable to black hole 
physics [8]. 

It is our aim to tidy up some dangling issues and fill technical details of the constrained 
system originated from the simplest version contained in the spirit of the initial proposal [I] . 
Based on those findings we will construct a sample of explicit cosmological vacuum models 
consistent with the new dynamics. Our final goal is to show in a cristal clear way that giving 
up full diffeomorphism invariance, in a very specific way, implies some bizarre and so far 
unexplored consequences at the most fundamental level. 

We will always compare our outputs with the deep infrared theory we deal with. General 
Relativity (GR) -with cosmological constant- [9], 

J^an = V7n(^''r + K,jK'^ -K^- 2A) , (1) 

where there is full 4-diffeomorphism invariance — > — e'^, with an arbitrary infinites- 
imal function of the coordinates as well as of the fields. The simplest setting one can find 
in the literature that mimics the partial breaking of the symmetry of time diffeomorphisms 
is just a small modification in the kinetic term of ([T|), see section 2. Notice however that 
even such a slight modification dramatically modifies the constraint's structure of the theory 
in phase space and makes it inconsistent in most classical settings of the initial conditions, 
with the sure exception of some symmetric, protected cases and of course by dimensionally 
reducing it to a one dimensional mechanical model; this shows how finely tuned is GR for it 
to be dynamically consistent. 

A general consequence of the modification of ([T]) - or the more general modifications 
discussed in the literature - is that, due to the reduction of diffeomorphism invariance, 
the foliation of the spacetime in equal-time surfaces is fixed and acquires a direct physical 
meaning: there is an absolute concept of simultaneity of events. In this sense the time in 
these modified theories of gravity is a Galilean time. 

As other consequences of this partial breaking of diffeomorphism invariance we may list 
some of our findings, see section H] 
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1. Due to the impossibility to perform time diffeomorphism with spatial dependences, one 
can not construct a Rindler causal horizon through each spacetime point and hence 
nor a Unruh radiation can be defined. 

2. Whereas the de Sitter solution in pure GR with cosmo logical constant accepts sev- 
eral types of time foliations (open, closed and fiat 3-slices) but the physics remains 
always the same, we find that in the modified theory of gravity, every type of foliation 
corresponds to a different physical vacuum. Hence the vacuum is degenerate. 

3. Some specific solutions of GR are easily seen to be solutions of the new dynamics as 
well, but one should be aware that this is not a coordinate independent statement. 
We will give examples of background configurations of GR that in some coordinate 
incarnation are indeed solutions of the minimally modified dynamics, whereas they are 
no longer so when written in other systems of coordinates. 

4. Contrary to some remarks in the literature, the implementation of the partial gauge 
fixing of the projectability conditions -i.e., that the lapse only depends on the time 
coordinate -, has no effect whatsoever on the Hamiltonian constraint, if one sticks to 
the rules of deriving the dynamics from the action principle. 

Due to a certain confusion in the literature, let us make an additional comment on the last 
remark above. One can interpret the projectability condition in two different ways. The 
natural way, which we advocate here, is to understand it as a partial gauge fixing for a 
theory defined by an action, with no preconceptions on the possible coordinate dependences 
of the fields. We may call it the soft projectability condition. Another view, the hard 
projectability condition, assumes a modification of the original action, in that it already 
contains the restriction that the lapse only depends on time. We believe that this second 
point of view, though legitimate, is less natural than the first one. We will be back to this 
issue in subsections 12.21 and 13.31 

Our analysis will be roughly split into two main venues, see section |3l 

i) In the generic case the constraint analysis can be performed to its end, and one obtains 
as a result a tertiary constraint and a partial determination of the Lagrange multiplier 
A° in the Dirac Hamiltonian. But one must be aware that other retrictions, coming form 
boundary conditions, may produce a collapse of the dynamics, in the sense for instance of 
the compulsory vanishing of the lapse [5]. 

ii) In non generic cases, we show that the typical situation is that of an indefinite chain 
of tertiary, quartiary, etc, constraints which seem to set notable limitations to the allowed 
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initial conditions. Even if we show some settings that allow for solutions, it is quite likely 
that there is inconsistency for a large class of initial conditions, much beyond the situation 
in GR. 

In addition of the mentioned bizarre consequences derived from partial diffeomorphism 
breaking and the increasing experimental evidence on the correctness of GR in the deep 
infra-red [TU] it is fair to mention that although in principle appealing, a weakness of the 
proposal is that there seems to be no clear mechanism guarantying that the low energy limit 
of this modified theory is going to restore full fledged GR. 

2 Set up: minimal deviation from General Relativity 

In order to make our point more clear we will follow the approach of [TT] which in some sense 
is based in some sort of effective theory construction. Such approach is theoretically robust 
and it is suitable from the phenomenological point of view. In doing so we will encounter 
two problems that we will refer below. The construction goes as follows: one introduces an 
anisotropic scaling at the ultraviolet fixed point from which one constructs a power-counting 
renormalizable theory organized in terms of a placeholder symbol k, which allows for a kind 
of weighted scaling dimensional analysis, 

[dx] = [K]-\ [dt] = [K]-\ [N] = [k]\ m = [K]\ [7d = M°- (2) 

At high energy this scaling modifies the kinetic part of ([1]) with the presence of a new 
parameter and simultaneously introduce a plethora of new relevant operators in the potential 
term. Most of these terms become irrelevant in the infrared. In most of applications one 
assumes that a naive power counting in derivatives holds and then one can discard the 
contribution of the tower of operators present in the potential in the infrared. Unfortunately 
their coefficients have a logarithmic running spoiling this behavior unless its ultraviolet 
value is unnaturally small [1], this is the first of the problem referred above. The second one, 
already stressed in [2], is that the inclusion of the Lorentz violation term in the kinetic part 
of the Lagrangian leads to an extra scalar mode for the graviton at all momenta and thus 
if one does not make corrections for these extra degrees of freedom one will never recover 
General Relativity (GR) in the infrared. A proposal has been made to amend such behavior 
[12] but some criticism on the results has been raised [13] . 

To avoid the first of these problems we will only focus in the possible deviations of 
GR near the infrared, thus in some sense we will give up issues concerning the ultraviolet 
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completion of the theory. With respect to the second problem we will not be concerned with 
it and explore only a minimal deviation from GR in the spirit of the original proposal [T]. 
Given the previous line of reasoning, the starting point see below, is the simplest of the 
settings discussed in [11] and corresponds to a set up where only relevant, dimension four 
operators in the infra-red are kept. The kinetic term will contain a slight difference with 
respect to GR while the potential term will remain the same. 

Let us remind of the 4-dimensional metric decomposition in terms of the 3-dimensional 
metric and the lapse and shift fields 
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(3) 



which implies y/—g = ^/^N. Spatial indices are raised and lowered with 'jij and its inverse 
7*-' respectively . The Lagrangian takes the form 

C = {^'^R + KijK'^ - XK^ - 2 a) , (4) 

(3) 

where R is the scalar curvature for the 3- metric jij, written henceforth simply as R, and 
Kij is the extrinsic curvature of the equal time surfaces, 

= ^ - V.iV, - V,iV.) , (5) 

with K been its trace, K = K^^'-fij. Notice that for the value A = 1 we recover the standard 
ADM Lagrangian ([T]) with all its virtues and properties. Thus the initial assumption in [T] 
is that as one goes to the infrared and recovers full diffeomorphism invariance the parameter 
A must run exactly to one. For A 7^ 1 the invariance is reduced to the sugbroup of foliation 
preserving diffeomorphisms, for which the time diffeomorphism variation dx^ =: — e*^ only 
depends on the time coordinate x^. From the phenomeno logical point of view the value of 
A is restricted to lie very near the GR value [H]. In this sense we will take A as a varying 
parameter that only in the deep infrared matches the GR value. For the sake of completeness 
we will explore also the region where A differs significanlty form 1. 

In order to obtain the Hamiltonian dynamics we define first the variables in phase space. 
Since the Lagrangian (jl]) does not depend on the time derivatives of the lapse and shift, we 
identify the momenta conjugate to the lapse =: and shift A^* as primary constraints 
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in phase space 

^ . (6) 
The Lagrangian definition of the momenta tt*-' conjugate to jij gives 

^(k'^ - XKfA . (7) 



Note that vr*-^ is a tensor density of weight 1 with respect to 3-diffeomorphisms. To prepare 
for the construction of the canonical Hamiltonian, we trade the canonical variables tt*-' for 
K'K The trace of ^ is (vr := n'^jij) 

7r = ^{l-3X)K. (8) 

At this point we notice that a special behavior takes place at A = I- In this case a new 
primary constraint appears, vr ~ 0. What happens is that, for this value of A, the gauge 
symmetry of the theory is enhanced: this is the anisotropic Weyl invariance. It is obvious 
that to a new gauge invariance there should be associated new first class constraints, as the 
eventual generators. We will leave aside momentarily this particular case and come back to 
it in section [5l 

Using ([7]) and ([8]) we can isolate K^^ in terms of the canonical variables, 

K^i = + tAtT^T^^) . (9) 

Vl 1 - 3A 

With all the previous inputs, we can rewrite the system in the Hamiltonian formalism. 
This will allow us to handle the stabilization - i.e. time preservation - of the constraints in 
phase space. The canonical Hamiltonian is defined by the spatial integration H = j d^l-L, 
with l-i = TT^^^ij — C. We obtain, up to boundary terms 

n = y/^N(^ -R + -f~\'^ni, + Y^^7^'^' + 2A) - 2N, Vjir'^ . (10) 

Defining (3 := 2(i-3\) ' Hamiltonian density can be cast as 

n = n'^''"' + p^n\ (11) 

which for (3 = reduces to the standard ADM formulation including the cosmological con- 
stant term. 
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2.1 The secondary constraints 

The secondary constraints are obtained under the requirement that the primary ones, P^, 
are preserved under the time evolution. By varying the action with respect to the lapse we 
obtain the Hamiltonian constraint 

•Ho = - i? + I-'tt'^'it.j + ^2 + 2a) , (12) 

while the variation with respect to the shift leads to the momentum constraints 

nj = -2\/in'^. (13) 

Notice that the only difference with respect to the ADM formulation lies solely in T-Iq. All 
in all we have 

'Ho = 'H7''+P^, H,=^r"^ (14) 

where {3 parametrizes the deviation with respect to GR. The Dirac Hamiltonian can be 
expressed as 

Hr>= ! di^{Nm^ + VP^) , (15) 



with the Lagrange multipliers A'^ being - as of now - arbitrary functions. 
2.2 The tertiary constraint 

Now we must look for tertiary constraints, as the consequence of the stabilization of the 
secondary ones ( !T4|) . "H^. This computation has been properly addressed in [3] within the 
Lagrangian framework in the tangent bundlJl]. Instead, we will reobtain the tertiary con- 
straint working within the canonical formalism. 

The stabilization of the secondary constraints takes the form of the requirement 



{n,, Hd} = !^n„ J d^'N'-H^i^')'^ 



(16) 



Using f|T5|) we can display the different terms that contribute to f|T6|) . Two relevant pieces of 
information help to calculate f lT6|) . They are: i) Notice first of all that due to the standard 



^Actually, from the Lagrangian point of view this wiU be a secondary constraint, because a n-ary constraint 
in phase space has a corresponding (n — l)-ary constraint - through the pullback operation - in the tangent 
bundle. Primary constraints in phase space correspond to identities in the tangent bundle. 



6 



behaviour of a scalar density under 3-diffeomorphisms one has 

7^0, f rfx'iV^Hi(x') j = diiN'-Ho) , (17) 



which vanishes in the primary and secondary constraints' surface, ii) On the other hand, 
the algebra of the generators of 3-diffeomorphisms closes in the standard way because they 
coincide with those of the ADM case. This in turn guaranties that the stabilization of the 
constraints Hi does not introduce new constraints. Using these two facts the only relevant 
term in the computation becomes 

tertiary constraint = |'Ho(x), J dx.' NHq^x.')^ . (18) 

The computation of (fT8|) is facilitated by the smearing of the constraint Hq by means of 
an arbitrary function t], which we take of compact support. So our task is reduced to the 
computation of 

j (ix 7] (tertiary c.) = rfx777{o(x), j (ix'iV?{o(x')| ■ (19) 

Using f lT4|) . we expand the rhs of f lT9|) as 

j (ix?7 (tertiary c. ) = C(/3°) + 0{/3) + Oif3^) . (20) 

The 0{/3^) term in the above expression is proportional to the momentum constraints and 
plays no role. The 0(/3^) term vanishes because no derivatives of the fields are involved. 
Thus the tertiary constraint stems only from the terms 0{f3). 

j dxr/ (tertiary c.) / dx'N'H:^^''"'^ ^ ^ d^N^, J dx' t]?/*^'''"' |^ .(21) 

Consider the first term in the r.h.s of (l2T]) . Getting rid momentarily of the smearing, we 

2 

realize that it corresponds to the time derivative of under the evolution provided by the 
ADM Hamiltonian in the gauge A^* = 

{^./dx'N^T'"}. (22) 

Thus we can compute it with standard formulas - see for instance [9l [15] - for the time 
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derivative of the components of the 3-metric and their conjugate momenta. We notice then 
the crucial fact that, when smearing with rj, all terms in f l22l) that have no derivatives of the 
fields will cancel against the second contribution to (12T|) . It turns out that 



( ) = irrelevant terms — 47r AA^ , (23) 

where AA^ stands for the Laplacian 7jjV*V-'iV. The relevant contribution in (l23l) arises 
entirely from the tt term. All in all, the contribution to (^T^ becomes 

j rfxT/ (tertiary c.) = -4/3 j d:si(r]7iAN - NnAr]^ =4/3 J d:s.r](2Vi7rd'N + NAi^j ,{24) 

where in the last equality we have used part integration and the fact that the function r] is 
of compact support. From the above expression the tertiary constraint is identified as 

2 ViTT (9W + A^Att ~ . (25) 

This equation, (1251) . has been first obtained in [16]. It has also been considered in [5], 
although not really interpreted as a constraint because the momenta canonically conjugate 
to A^^ are eliminated in [3] and the lapse and shift variables take over the role of Lagrange 
multipliers. We prefer for now to keep all the variables of the formalism and the Hamiltonian 

m- 

At this point it is worth noticing that from the perspective of the hard projectability 
condition, there is no tertiary constraint at all, because in such case one must stabilize 
j (ixHo, which is like redoing with r] = 1; then the computation in the first equality 
of (l2l|) gives a vanishing result: rjirAN — NnAi] = 0. In conclusion: there is no tertiary 
constraint for the hard projectability condition implemented in the Lagrangian (jlj). 



3 Stabilization of the tertiary constraint 

The stabilization of fl2^ under the time evolution will produce the appearance of N, which, 
according to the dynamics, is the multiplier A" in the Dirac Hamiltonian (fT5|) . Thus the 
stabilization of (125!) has the form 

Att A° + 2Vi7r a'A° + (terms with no A°) = , (26) 
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where we have taken the equahty to zero as an ordinary equahty because it is legitimate to 
make the determination of the multiphers A'^ on shell. Note that there are several possible 
choices in the space of field configurations that lead to the fulfillment of fl25|) and fl26|) . 
For instance, if we consider a configuration satisfying VjTr = 0, then fl25|) holds and fl26|) is 
independent of the value of A°, but not void. 

In this section we will consider different consistent ways for which the tertiary constraint 
is enforced by restricting the field configurations in a way that guaranties ( p5|) . 

3.1 The generic case: Att ^ 

As a first trial let us consider a generic case with Att 7^ 0. Notice that if i) the initial 
conditions for the fields at, say, t = 0, satisfy all the constraints, including fl25|) . and if 
ii) the multiplier A° in the Dirac Hamiltonian f lT5]) is chosen so that it satisfies f l2B]) . then 
it is guarateed that all constraints will be satisfied at any time. Which means that the 
constraint analysis is finished. Notice though that the analysis made here is based on local 
requirements. What may remain to be explored is the adequacy of these results for certain 
boundary conditions. Some concerns in this respect will be drawn below. 

In this generic case one can obtain, at least formally, the partial determination of the 
multiplier A° as follows: First factorize out Att from fl26|) and obtain an expression of the 
form 

AO = yA° + t/ , (27) 

with V the differential operator V = —2 / di and U the remaining quantity. The solution 

Att 

to (!27|) will be the sum of a particular solution plus an arbitrary solution of the associated 
homogeneous equation. As long as the vector field V is diferent form zero - which must 
be, because we are in a generic configuration and so Vtt 7^ -, one can locally change to a 
spatial coordinates system such that V = dz- Then the general solution of fl271) is 

A°(a;, y,z,t) = -j dz'e^~^'U{{x, y, z' , t)) + e^f{x, y, t) . (28) 

As already mentioned, our determination of A*^ is only of local validity and one may 
ask on the global behaviour of the function A*^, which is an issue not addressed here and 
that essentially depend on the specifics of the boundary conditions that are imposed. This 
behaviour will affect that of A^, since AT = A°. Within the very same generic case and 
assuming an asymptotically flat space the behaviour of A^ has been studied directly from 
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the analysis of (125|1 in [5]. The outcome of the analysis is that to prevent the lapse function 
from been divergent at the boundary, the only acceptable solution is = 0. 

The cases to be considered in the following subsections are non generic, that is, they 
satisfy Avr = 0. 



3.2 Preserving the first class condition. 

Let us consider a second trial set for the fulfillment of the tertiarty constraint (125|) . Since 
the consistency of the whole picture requires that 7^ 0, otherwise the Lagrangian vanishes 
and the dynamics disappears, we can take the tertiary constraint in the form 

$ := 2Vi^ + Avr ~ , (29) 

which is more suited for our present purposes. Note that except for some special circum- 
stances, this constraint will make the primary constraint Pq := P second class, 

{P(x), $(y)} = ^(y) (^(y)<5(x - y) - 2Vr5(x - y)) ^ , (30) 



even in the case when the projectability condition holdcl unless W vanishes on shell. This 
last situation amounts to impose the new set of constraints 

Vivr ~ , (31) 

which already imply (!29|) . As a matter of fact, there is another reason for analyzing this 
case: if we were to require that (125|) introduces no restriction upon the lapse, then we should 
impose (1ST]) . Since the partial spatial derivatives commute with the dynamics, it is more 
convenient to use the form 

(32) 



in order to explore the stabilization of ( 1311) . After using the Hamiltonian constraint to 
eliminate terms with n'^^nki, the new constraints originated from the stabilization of fl32|) 
become 



N [2R 6A - 2AA^ 

1-3A 7 



0, (33) 



^Notice however that being second class in field theory is more tricky than in the mechanical case. Indeed, 
if we were in mechanics, a second class condition for P will undoubtedly determine the multiplier A°. As we 
have just seen, this is not the case in field theory. 



10 



or, since di{ — ) vanishes on shell by virtue of 

7 



^ - TTT^ ^ 3A + - ~ , (34) 



V 2(1-3A) 7 ) N 

thus showing that the stabihzation of (13T|) has introduce a new constraint flM|) . Note that 
the initial goal to prevent the constraint Pq to become second class faces new challenges. It 
seems that the most simple setting aimed at this contains as an ingredient the projectability 
condition, i.e. the lapse only depends on the time coordinate, together with the requirement 
that the 3-surfaces - labeled by the time coordinate - of the foliation of the spacetime must 
be surfaces of constant curvature, that is, 

~ , diRc^O. (35) 

But we are not over yet, because, again, stability must be required to these constraints. In 
this new analyisis we must consider i) The stabilization of the projectability condition in (155]) 
results in the multiplier A" satisfying diX^ = on shell, which is just fine because the only 
time-diffeomorphisms permitted are those of the type — )■ — (5ge'^(x*^). ii) As regards 
the stability of diR ~ 0, notice that since i? is a 3-scalar, we know that {R, J N'^'Hi} = 
N^diR, which already vanishes on shell by virtue of (1551) . Thus we only need to compute 
{diR, J NTio}. This results in the new set of constraints 

a.(^^)^0. (36) 

Other constraints will follow from the stabilization of ( |36|) . We will not pursue this route any 
further, but make a general comment. There seems to be no obstacle for the stabilization 
mechanism originated from ( 1321) to go on indefinitely, giving at every stage new constraints 
of the form of the gradient of a scalar. Notice that this scalar will involve more derivatives 
the more we advance in the process. 

3.3 The projectability condition 

The third and last scenario we will focus on for the fuUfilment of the tertiary constraint 
([25]) assumes that the projectability condition is satisfied. In this long are we keep 
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7^ 0, a new constraint is compulsory, namely, 



Att ~ . (37) 

The stabilization of the projectability condition goes along the very same lines as in the 
previous case, subsection 13.21 As regards ( 1371) . the easiest way to satisfy it is to asume the 
stronger condition f l32|) . in which case we will end up with the scenario discussed in 13.21 
Going back to the stabilization of f l37|) . and following the same argument as with fl32|) . we 
will stabilize 

v/7A(-^)^0. (38) 



Since A(^) is a 3-scalar, we already know that {^A(^), / Nmi} = di(^N^A{^)^ ~ 0, 
and we must only compute {^y^A(-^), J NT-Lq}. Keeping 7^ the final result is 

(v^A(^))" ^ 2N^^AR - d,{n^^d,{^))) , (39) 



/7V V- ' - -^7^ 

and thus the new constraint, byproduct of the stabilization of (138|) . is 

^AR-d,{n'^d,{^))c^O. (40) 

As in the previous case the stabilization mechanism goes on indefinitely, and an infinite tower 
of new constraints, with more and more derivatives involved, appears as we advance in the 
algorithmic procedure. 

In some versions of the original model one implements the projectability condition from 
the very begining. It is claimed that within this proviso the Hamiltonian constraint is not a 
local equation satisfied at each spatial point but an equation integrated over a whole space. 
This observation is misleading in one respect: when one implements a gauge fixing (which is 
the spirit of the projectability conditions) within the action itself, one may lose information 
on the constraint structure of the theory [T7j, which must be restored by hand, and the 
information lost in this case is precisely the local structure of the Hamiltonian constraint. 
Of course one can always change the theory and consider that the projectability condition is 
a direct ingredient of the action principle, in which case the nonlocal Hamiltonian constraint 
appears, but this runs againts the spirit of writing a Lagrangian with no preconceptions as to 
whether its dynamical consequences could be. Sticking to an action principle with no further 
preconceptions, one can introduce the projectability conditions as a partial gauge fixing, as 
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done above, but it does not have consequences regarding the localhty of the Hamiltonian 
constraint. 

Before closing this section, let us make three relevant remarks concerning the previous 
two subsections, 13.21 and 13.31 First of all, both outcomes are similar to that encountered by 
|18] . where it is argued that an infinite set of constraints appears in the theory, but whereas 
in [18] the focus is in the potential term and A is kept to its GR value, our analysis is for 
A 7^ 1. In spite of the existence of an infinite chain of constraints, particular configurations 
compatible with the full set of constraints exist, as shown in subsections 14. II and 14.21 below. 

The second remark, concerning a possible way out of this infinite string of constraints, 
is the extreme, but consistent setting, of requiring that the variables jij , vr*-^ ,N,N^ be 
only time dependent. This outcome is nothing but the dimensional reduction of to a 
mechanical model with only one coordinate - the time. The reduction is anticipated to be 
consistent because it is made under the abelian group generated by the Killing vectors di. 
The only gauge freedom left is that of time reparametrizations. Though consistent, this 
extremely restrictive framework seems of very limited interest, at least for what regards its 
classical field theory content. Obviously, for /3 — we do not recover GR. 

Finally, the last remark concerns the issue of consistency. We expect that the addition 
of higher ordre terms to the minimal setting discussed here, (jlj), does not change any of our 
conclusions. In fact, it will just add further complications to the constraint analysis, which 
will appear even more restrictive. 

4 Some well-known solutions satisfying the constraints 

To complete our findings we give some explicit examples of solutions. As such, they fulfill 
the condition (125|) and as a consequence any possible chain of constraints derived form its 
stabilization. We choose two different set ups, subsections 14.11 and 14. 2t the first one does not 
imply the projectability condition whereas the last one does. 
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4.1 Black hole solutions 

An obvious setting which guaranties f l25p is to require the vanishing of TTjj which in turns 
imphes, recalhng the Hamiltonian constraint f lT2|) . 

~ 2A . (41) 

Notice that tt*-' ~ directly imphes K^^ ~ 0, and hence the parameter A in (jl]) and in 
the equations of motion (eom) (IA.1|) . flA.2|) plays no role at all, as it is directly seen in the 
Hamiltonian (fTUj) . 

4.1.1 Schwarzswild black hole 

By direct inspection of (|5]) we can conclude that K^^ ^ is satisfied by any stationary 
solution with vanishing shift A^*. In such case we need to explore the stabilization of the new 
constraints 

A^* ~ , TT^J' ~ . (42) 

Preservation in time of the vanishing of the shift trivially determines the arbitrary functions 
A* = in the Dirac Hamiltonian, while the stabilization of tt*-' ~ gives the new constraint 

NR., - {V^VjN - ^ijAN) ^ . (43) 

The trace of (143|) implies, using (14T|) . the following constraint on 

AA^ + A^A ~ , (44) 

which is no longer compatible with the projectability condition, as long as we keep N ^ 
and A 7^ 0. Finally, inserting ( I44p in ( I43p . the new constraint is 

NRij - ViVjN - ANji^ ^ , (45) 

which already implies fl44p . if one takes into account 041 p . The stabilization of fl45p yields 
a partial determination of the arbitrary function A° in the Dirac Hamiltonian but (145 p can 
also be seen as a partial determination of the lapse. 

Intepreting equation as an equation for the lapse, Rij must satisfy some projectabil- 
ity conditions. Let us find them. Computing the divergence of and using the result 
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VjViVN = diAN + RijV^N one obtains the equation 

NV^Rij - di{AN + NA) ~ . (46) 

The last piece vanishes using the constraint itself, since ( 1^ and (145|) imply (jH]). Factor- 
izing the lapse, which we assume is a non-vanishing function, we arrive at the integrability 
condition for fH5|) 

V^R,j - . (47) 

Using the contracted Bianchi identity, this condition is equivalent to Vji? ^ 0, which is 
satisfied because of 0411) . Thus the integrability condition is satisfied, showing the consistency 
of our procedure. 

As an application of this setting, consider the static, spherically symmetric ansatz 

ds^ = -A{r)df + B{r)dr^ + r^dQl (48) 

This ansatz has been already considered in [5] but under the detailed balance condition. Eq. 
f l48p is in fact a consistent reduction, which sets the shift vectors to zero as well as eliminates 
other metric components. The solution fulfilling the eom and the constraints fj^Tl) and fj45|) 
boils down to just Schwarzschild, either de Sitter or Anti de Sitter 

2 A 

A{r) = B-\r) = l-^ + -. (49) 

4.1.2 Kerr blak hole 

One may wonder whether any solution in GR admits a continuation to a solution in the mod- 
ified theory (jl]). It is clear that solutions to (jl]), which must fulfill the constraint 025p . yield 
in the limit A — 1 solutions of GR. But precisely, because of the existence of the additional 
constraint ( 125|) . the opposite is not true. Candidate GR solutions that, in principle, may 
fail with the fulfillment of (125|) are those with non- vanishing momenta ir^^ . If the solution 
is stationary this means that the shift should be different form zero. A natural candidate is 
therefore a rotating black hole, and a possible choice inside GR is the Kerr black hole - we 
set A = for simplicity. This solution will have a continuation to a solution of (jlj) as long as 
f l25|) is satisfied. Let us check whether this is the case. Writting the metric in the dragging 
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frame, the non-vanishing components of the 3-metric are 



n2 



P 2 
9rr — -r , gee — p , 9^ip 



((r^ + a^) + ^^^^^^)sin^^, (50) 



with Tg the Schwarzschild radius and 

J 2 2,22/) A 2 ,2 /rn \ 

a=— , p=r+a! cos fc* , A = r — rr^ + a . (51) 
Using ( l50l) we can obtain the 3-dimensional shift and lapse fields 

/V2 - _i , '^^^ rV.^g^sin^^ 

p p"^(p^(r^ + a^) + rr^ a sm u) 

AT = -(0,0 



+ a2)p2 + rr^a^ sin ^ 



One can check from (!50|) and (152!) that vr*-' 7^ 0, but its trace vanishes, tt = 0. This means 
that (!25|) is still satisfied. Since one can read from f lA.l|A.2p that the difference between the 
eom of GR and those of (jl]) lies solely in terms that are proportional to vr, we can immediately 
conclude that the Kerr solution fl50|) is also solution of (jl]) for any value of A. 



4.1.3 The Galilean nature of time for A 7^ 1 

In view of the previous two examples one may be tempted to erroneously conclude that any 
solution to GR with sufficient symmetry can always be deformed in the A-parameter space 
to fulfill the new dynamics defined by (jl]). In fact, for the theories defined with (jlj) with 
A 7^ 1, the restriction to foliation preserving diffeomorphisms has the consequence that there 
is a preferred frame with Galilean time, thus restoring and absolute concept of simultaneity. 
The assertion that Schwarzschild, or Kerr, are solutions of (jlj) must be qualified: they are 
solutions in the coordinate frames that only differ in a foliation preserving diffeomorphism 
from the expressions given in the last subsections for such solutions. But if these solutions, 
as solutions of GR, are presented in other frames, they are no longer solutions of (jl]). Let us 
give the example of Schwarzschild in - ingoing - Painleve-GuUstrand coordinates. 



ds"^ = -dt^ + {dr + v^l - A{r)dtf + r'^d^l , (53) 
where the 3-metric is simply flat, the lapse is trivial and the shift has the only component 



A^^ = y/l - A{r). One can immediately check that Att 7^ 0. This result, in view of A^ = 1, 
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shows that the constraint (125|) is not fulfilled and therefore (153|) is not a solution of @. 
This result above is not in contradiction with being a solution of (jlj, because the two 
types of coordinatization used in psjl and fl53l) are not connected by a foliation preserv- 
ing diffeomorphism, and therefore the backgrounds described with these coordinatizations, 
though they lead to a single physics when we go to the GR limit, represent different physical 
settings for the dynamics given by (jl]). Later we shall complement these findings in terms of 
Rindler coordinates and examine its immediate consequences for thermodynamics, For the 
time being let us make the remark that one can associate a Hawking radiation in both sys- 
tems of coordinates, and thus this emission is insensitive to the partial time diffeomorphism 
breaking. 

4.2 Vacuum cosmologies, generalities 

A less restrictive non generic setting, still fulfiling ( 125|) . is constructed by demanding 

n'^ = b{t),/^Y^ , (54) 

in the gauge A^* = . We will obtain time dependent backgrounds that will lead to a number 
of interesting consequences for well known cosmologies. As a matter of notation, and since 
we are directly looking for solutions, we write henceforth ordinary equalities when dealing 
with the constraints. 

First of all notice that the geometry of the 3-surfaces implies that b in fl5^ is a 3-scalar, 
because tt*-' is a tensor density; thus the assertion that b depends only on the time coordinate 
is a covariant statement with respect to 3-diffeomorphisms. And it is also so with respect to 
the time foliation preserving diffeomorphisms. 

Inserting (154|) in the eom (lA.ip and solving for 7 one obtains 

/ 2A /■* \ 

7ij(x, t) = exp I _ b{T)dr \ 7ij(x, 0) , (55) 

that is, the evolution of the 3-metric in terms of expansion factor, exp(Y^^ j^b{T)dT), and 
the initial data. Notice that this expansion is rigid in the sense that all points on the surfaces 
evolve with the same factor. 3-dimensional, distances between them are only affected by a 
global time dependent factor. 

In the sequel we will obtain the dynamical evolution of b{t) for different cases. Let us 
stress nevertheless that this is not the most efficient way of constructing time-dependent 
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metric spaces which are solutions of (jlj) should one be interested in including matter. 
Notice, first of all, that with the ansatz f l54p the Hamiltonian constraint (125 p becomes 

R=-J-^b\t) + 2A, (56) 

which enforces R to depend only on time. From the trace of (1541) we infer that the tertiary 
constraint (125|) is satisfied because VjTr = 0. Note that our ansatz complies with the assump- 



tions of subsection 13. 2[ including the projectability condition of the lapse, to be obtained 
below. 

There are two ways to compute the time derivative of tt*-' : either we use (IA.2p or the rhs 
of the ansatz (154|) . Using the former we obtain 



( - + li?7*^' - i' - A7*^) , (57) 



2(1 - 3A) 

while the latter leads to 



= J^Ml'') = {b+ Y^b') . (58) 



Equating both expressions, taking the trace, and using ( l56l) . we obtain our fundamental 
equation 

We will explore the solutions to fl59|) later on. Note that fl59|l enforces the projectability 
condition on the lapse and it also indicates that n^^ is proportional to the metric Applying 
this result to (IFTp . and using fISBl) and fISIJl) . we obtain 

We infer from equation ( 160|) that the surfaces of the foliation are Einstein spaces. The 
36 

3-curvature, ^ ~ ~ opposite sign to b. 

According to the signs of the cosmological constant A and the coefficient 1 — 3A we will 
have four different cases that we consider below. In all the following set ups, A^ is gauge 
fixed to 

|1-3A| , , 

N=^-^-^, (61) 
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and a > is defined as 

\ ■■= . (62) 

4.2.1 de Sitter-like space: 1 - 3A < 0, A > 

This case contains GR - with positive cosmological cosntant - as a particular case and seems 
to be phenomenologically the most plausible scenario [HI [19]. Then the equation (159|) 
becomes 

(63) 

Its solutions are separated in three regimes, accordingly with the initial conditions, given by 
2 

the sign of \b\. 

a 

Closed slicing. If the initial data fulfills the inequality 

|b|<-, (64) 
a 



the solution to ( 163|) is given by 



h{t) = --tanh ( - ) , (65) 



a \a 



up to a rigid time translation. For this solution b is negative. The simplest way to 
realize the 3-dimensional Einstein space with positive curvature is the 3-sphere which 
has Rij = 2jij. Using ( l55ll . the 4- metric becomes 



ds^ = -NHt" + Na' cosh^ [^]dnl. (66) 



Open slicing. If the initial condition is such that 



|6|>-, (67) 



a 



the solution to ( 163|) reads, up to a rigid time translation 



h{t) = --coth ( - ) , (68) 



a \a 



and now b is positive. The simplest way to realize the 3-dimensional Einstein space 
with negative curvature is the 3-hyperboloid, which has Rij = —2jij. With a procedure 
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similar to the previous case, using (155|) . the 4- metric becomes 

ds'^ = -N'^dt^ + Na^ sinh^ dH^ . (69) 

This case admits the hmit a — ?■ oo, where fl69p becomes 

_ _N^clf + N fdHl , (70) 

which, for A = 1, reduces to Milne universe, i.e. FRW solution in absence of matter. 
Flat slicing. As a last case we encounter the critical situation 

\h\ = -. (71) 
a 

This case has the trivial constant solutions 

h{t) = ±^ . (72) 

The 3-curvature vanishes and we realize it with a spatial flat euclidean space. The 
4-metrics are 

ds^ = -N^df + exp (^±^^) (dx^ + dy^ + dz^) . (73) 

The solutions (l66l) . ( l69|) are only solutions of GR for A = 1. Instead, ( l73l) is a solution of 
GR for any A if one rescales conveniently the cosmological constant. All the above solutions 
(!66|) . (!69l) . (!73|) . boil down to the same GR solution for A = 1: de Sitter space. For A 7^ 1 
the gauge group is reduced to foliation preserving diffeomorphisms, and as a consquence 
the three cases, open, closed and fiat slicing, are indeed three different vacuums. This is 
vacuum degeneracy. Taken this result at face value, if nature had initially broken the 4- 
diffeomorphism group in this region of the (A, 1 — 3A) plane, it would also had to undergo an 
spontaneous symmetry breaking [22] in order to conform with the present-day observationally 
allowed values of A and A. This mechanism of spontaneous symmetry breaking is also present 
in other models of Lorentz violation with an unconventional kinetic term [23]. Althoug this 
picture is appealing, one should bear in mind that the analysis made here is classical, without 
including matter, and susceptible to quantum corrections. 
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4.2.2 Anti de Sitter-like space: 1 - 3A < 0, A < 



The equation (139]) becomes 




(74) 



which has as solution 




(75) 



up to a rigid time translation. Now b is positive and the 4-metric becomes, 



ds"^ = -N^dt^ + Na^ cos^ - dHl . 




(76) 



For A = 1 this is just a partial covering of Anti de Sitter spacetime with coordinate singu- 
larities located at cos(-) = 0. 

Hitherto we have discussed cases that presume a mild modification with respect to GR 
in the sense that the value of |1 — 3A| in the (A, 1 — 3A) plane lies in the same, lower half 
part where GR resides. In that respect the remaining two cases we will discuss can support 
deviations far from GR because they lie in the upper half plane. 

It is worth noticing that if b{t) solves ( l59l) for some signs of (1 — 3A) and A, the change 
b{t) — —b(t) produces a solution of (159|) with the opposite signs, whereas the expansion 
parameter in (!55|) does not change because of compensating factor signs. This fact allows us 
to obtain immeditately the new solutions below from the previous ones above. 

4.2.3 1 - 3A > 0, A > 

This case is parallel to that in I4.2.2l but with the change b{t) — —b{t). The 4-metric becomes 




(77) 
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4.2.4 1 -3A > 0, A < 

Now this case is just that of l4.2.1l with the change b{t) —b(t). All previous considerations 
hold and we end up with the 3-metrics. 

ds^ = -N^dt^ + Na^ cosh^ dH^ . (78) 

ds^ = -N^dt^ + N sinh^ (^^^dnl 

ds^ = -N'^df + exp (±-t] {dx^ + dy"^ + dz^) . 



The considerations in the above section concerning the space of vacuum solutions can be 
summarized in figure [H 

A 



a 

-N'df + N a' sinli' { - ] d^ 



-N'dr + exp ( ±-1 ] {dx' + dy' + dz' 



N^dt^ + Na' cos" 



a 



dHi 



N'^df + Na' cos' 



a 



dQi 



►A 



ds' = 
ds' = 
ds' = 



-N'dt' + No' cosh' [-]dQ 



-N'dt' + N a' sinh' I - ] dW 



-N'dt' + exp ( ±-t ] (dx' + dy' + dz' 



Figure 1: Pictorial representation of the vacuum solutions in the (A, 1 — 3A) plane. The red 
axis represents the critical case A = -. stands for A^ = and — := -A^|A| 



Notice that while the transition between left f-)- right half planes as a function of A is 
smooth this is not the case for the top^botton half planes as a function of 1 — 3A. Notice 
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also the existence of a "quasi-inversion" symmetry by which we send 1 — 3A —(1 — 3 A) 
and A — —A and the spherical and hyperbolical 3-dimensional spaces are exchanged. 



4.2.5 The Galilean nature of time, revisited 

One of the disturbing outcomes of the model (jl]) is Lorentz violation, expressed by the fact 
that only time foliation preserving diffeomorphism are permitted. As mentioned earlier l¥. 1.31 
this is connected with the existence of a preferred time frame. Even if it has been advocated 
in cosmology that there is no reason to refuse the existence of that frame and furthermore 
that there exists such a natural, preferred frame defined by the cosmic microwave background 
[20] it is also true that, at least in the case of GR, tight phenomenological constraints rule 
out the existence of Lorentz violation operators of dimension < 6 [2T]. Thus it seems 
that consistency with phenomenological results does not leave too much room for Lorentz 
violation, unless it comes in a very exotic manner. Anyhow we will explore some of the 
possible consequences of Lorentz violation at most fundamental level. Our starting point is 
a solution of GR, the Milne universe, f l76p which covers the complementary wedge of Rindler 
spacetime in Minkowski spacetime. It is well known that, for t > 0, the following change of 
coordinates 



which is plainly just half of Minkowski spacetime, T > 0. We have gone with such detail, in 
these trivial matters, because of the points we want to make. Notice that for A 7^ 1, we can no 
longer practice arbitrary reparametrizations, but only foliation preserving diffeomorphisms. 
In particular, a change of coordinates of the type fl7^ is no longer permitted. This difference 
with GR is crucial in two essential points with a common origin: 

Cosmological time. In the GR case, the hyperboloid nature of the 3-surfaces of the foli- 
ation in fl69|) was just an artifact of the choice of coordinates for de Sitter spacetime. 
But in theories with A 7^ 1 we are stuck with such a foliation with hyperboloids, which 
are equal time surfaces, and the parameter t is already the physical time, except per- 
haps for a reparametrization involving only this time parameter itself. In other words, 
what for GR was just a coordinate singularity for t = in ( 169|) has become for A 7^ 1 
theories a true background singularity, a big bang for the solution with t > or a big 



R = t sinh ip , T = t cosh -0 



(79) 



makes this solution to be written as 




(80) 
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crunch for that with t < 0. 

In this spirit the three types of de Sitter-hke spacetime fohations can be interpreted 
as: Closed slicing ( l66l) describes a bounce. Open slicing ( l69l) describes a big bang or a 
big crunch with the singularity at finite time. The case of Flat slicing (173|) describes 
an infinite expansion or contraction with the singularity at infinity. 

Thermodynamics. Unlike Hawking radiation effect, Unruh effect is tantamount to Lorentz 
symmetry: in the absence of the latter the former does not exist [21]. The algebraic 
proof, in the framework of axiomatic field theories, states that the Minkowski vacuum 
restricted to a Rindler wedge is a thermal state with respect to the boost parameter 
[25]. It is obvious that in the absence of boost transformations of the Lorentzian type, 
we are not allowed anymore to define a thermal state i.e. one can not even find a 
horizon with a surface gravity from which to define the temperature through Tolman's 
law. 

Notice that besides these discrepancies with GR, the structure of the time diffeomorphism 
breaking allows the preservation of the continuos self similarity property on the solutions of 
the modified dynamics. 

5 The critical case A = | 

As anticipated in section [2l the critical case needs special attention because of the appearance 
of a new primary constraint, tt ^ 0. This means in particular that all solutions of (jlj), 
including the GR case A = 1, that happen to satisfy tt = 0, will be solutions for this 
critical case. And vice-versa: any solution of the critical case is a solution of (j4]) for any A. 
The canonical Hamiltonian may be taken as the same as in GR, because the A dependent 
term in ( ITOl) vanishes for tt ~ 0. In fact this term is quadratic in the constraint. The 
canonical Hamiltonian has always the ambiguity of the addition of terms linear in the primary 
constraints. The choice of the GR form for f|T2l) is the most convenient because of the 
closedness of the algebra of the Hamiltonian and momentum constraints in this case. This 
choice is legitimate because what is required to the canonical Hamiltonian is that {FL*p)q — 
FL*^ = C, where FL* is the puUback operation of the Legendre map FL from the tangent 
bundle to the cotangent bundle, and this condition is satisfied in the critical case for any 
value of A in ( ITOj) . 
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With this choice of the canonical Hamiltonian, the Dirac Hamiltonian is now 




{ADM) 



(81) 



Whereas the stabihzation of the primary constraints ^ yields the usual Hamiltonian 
and momentum constraints of GR, that of the new primary constraint vr ~ produces the 
new secondary constraint 



that translates in a partial determination of the lapse. Finally, the stabilization of the 
Hamiltonian constraint yields a partial determination of the arbitrary function ^, 



One can stabilize the constraint (182 p . giving a partial determination of A°. Alternatively one 
could have eliminated the momenta from the formalism and take the lapse and shift as 
the arbitrary functions of the dynamics for the 3-metric. 
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A Appendix 

In this short appendix we display for the sake of completeness the eom for the minimal 
dynamics dl]). Let us we stress that dotted quantities in subsections 13.21 and 13.31 stand for 
the evolution under this minimally modified dynamics and they should not be confused with 
dotted quantities used in subsection 12. 2[ which meant evolution under the ADM dynamics 
in the gauge A^* = 0. 

The eom for the dynamics (jl]) are just a minimal modification of that of GR [TJ] 



N{R - 3A) - AAT ~ 



(82) 



^(/? - 3A) - = . 



(83) 



_ 



1-3A 



A 



(A.l) 



25 



References 

[1] p. Horava, "Quantum Gravity at a Lifshitz Point," Phys. Rev. D 79 (2009) 084008 
[arXiv:0901.3775 [hep-th]]. 

[2] C. Charmousis, G. Niz, A. Padilla and P. M. SafRn, "Strong coupling in Horava 
gravity," JHEP 0908 (2009) 070 [arXiv:0905.2579 [hep-th]]. 

[3] D. Bias, O. Pujolas and S. Sibiryakov, "On the Extra Mode and Inconsistency of 
Horava Gravity," JHEP 0910 (2009) 029 [arXiv:0906.3046 [hep-th]]. 

[4] R. lengo, J. G. Russo and M. Serone, "Renormalization group in Lifshitz-type 
theories," JHEP 0911 (2009) 020 [arXiv:0906.3477 [hep-th]]. 

[5] M. Henneaux, A. Kleinschmidt and G. L. Gomez, "A dynamical inconsistency of 
Horava gravity," arXiv:0912.0399 [hep-th]. 

[6] G. Calcagni, "Cosmology of the Lifshitz universe," JHEP 0909 (2009) 112 
[arXiv:0904.0829 [hep-th]]. 

[7] R. Brandenberger, "Matter Bounce in Horava-Lifshitz Cosmology," Phys. Rev. D 80 
(2009) 043516 [arXiv:0904.2835 [hep-th]]. 

[8] H. Lu, J. Mei and C. N. Pope, "Solutions to Horava Gravity," Phys. Rev. Lett. 103 
(2009) 091301 [arXiv:0904.1595 [hep-th]]. 

[9] R. L. Arnowitt, S. Deser and C. W. Misner, "The dynamics of general relativity," in 
Gravitation: an introduction to current research, Louis Witten ed. (Wilew 1962), 
1962. arXiv:gr-qc/0405109. 

[10] R. Reyes, et ai, "Confirmation of general relativity on large scales from weak lensing 
and galaxy velocities," Nature 464 (2010) 256. 

[11] T. P. Sotiriou, M. Visser and S. Weinfurtner, "Phenomenologically viable 
Lorentz-violating quantum gravity," Phys. Rev. Lett. 102 (2009) 251601 
[arXiv:0904.4464 [hep-th]]. 



26 



[12] D. Bias, O. Pujolas and S. Sibiryakov, "A healthy extension of Horava gravity," 
arXiv:0909.3525 [hep-th]. 

[13] S. Mukohyama, "Caustic avoidance in Horava-Lifshitz gravity," JCAP 0909 (2009) 
005 [arXiv:0906.5069 [hep-th]]. 

[14] S. Dutta and E. N. Saridakis, "Overall observational constraints on the running 
parameter A of Hofava-Lifshitz gravity," arXiv: 1002.3373 [hep-th]. 

[15] See for instance Appendix E in R. M. Wald, "General Relativity," Chicago Univ. Pr. ( 
1984) 49 Ip 

[16] A. A. Kocharyan, "Is nonrelativistic gravity possible?," Phys. Rev. D 80 (2009) 
024026 [arXiv:0905.4204 [hep-th]]. 

[17] J. M. Pons, "Plugging the gauge fixing into the Lagrangian," Int. J. Mod. Phys. A 11 
(1996) 975 [arXiv:hep-th/9510044]. 

[18] M. Li and Y. Pang, "A Trouble with Hofava-Lifshitz Gravity," JHEP 0908 (2009) 015 
[arXiv:0905.2751 [hep-th]]. 

[19] F. S. N. Lobo, T. Harko and Z. Kovacs, "Solar System tests of Horava-Lifshitz black 
holes," arXiv:1001.3517 [gr-qc]. 

[20] S. Kanno and J. Soda, "Lorentz violating inflation," Phys. Rev. D 74 (2006) 063505 
[arXiv:hep-th/0604192]. 

[21] D. Mattingly, "Modern tests of Lorentz invariance," Living Rev. Rel. 8 (2005) 5 
[arXiv:gr-qc/0502097]. 

[22] R. N. Sen and C. Weil, "Degenerate vacuum formalism for spontaneous symmetry 
breakdown," J. Phys. A 4 (1971) 632. 

[23] See for instance N. Arkani-Hamed, H. C. Cheng, M. A. Luty, S. Mukohyama and 
T. Wiseman, "Dynamics of Gravity in a Higgs Phase," JHEP 0701 (2007) 036 
[arXiv:hep-ph/0507120]. 

[24] D. Campo and N. Obadia, "Why does the Unruh effect rely on Lorentz invariance, 
while Hawking radiation does not ?," arXiv:1003.0112 [gr-qc]. 



27 



J. J. Bisognano and E. H. Wichmann, "On The Duality Condition For A Hermitian 
Scalar Field," J. Math. Phys. 16 (1975) 985. "On The Duality Condition For 
Quantum Fields," J. Math. Phys. 17 (1976) 303. 



28 



